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Abstract 

In 2004 J. Fejoz [7], completing investigations of M. Herman's [9], gave a complete proof 
of "Arnold's Theorem" [1] on the planetary many-body problem, establishing, in particular, 
the existence of a positive measure set of smooth (C°°) Lagrangian invariant tori for the 
planetary many-body problem. Here, using RuBmann's 2001 KAM theory [16], we prove 
the above result in the real-analytic class. 
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1 Introduction 



The planetary many-body problem consists in studying the evolution of (1 + n) 
bodies (point masses), subject only to the mutual gravitational attraction, in the case 
one of bodies (the "Sun") has mass mo considerably larger than the masses rrii of 
the remaining n bodies (the "planets"). The Newtonian evolution equations for such 
problem (in suitable units) are given by 

m i <ij = m i m k — %r , i = 0, 1, n , (1) 



* Acknowledgments. We are indebted with Jacques Fejoz for many enlightening discussions. This 
work was partially supported by the Italian MIUR project "Metodi variazionali e equazioni differenziali 
nonlineari". 
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where qj = qj(t) € M 3 denotes the position at time t of the j th body, '| • |' denotes 
the euclidean norm and ' ' ' denotes time derivative. 

In [1, Chapter III, p. 125], V.I. Arnold made the following 1 : 

Arnold's statement: "In the n-body problem there exists a set of initial conditions 
having a positive Lebesgue measure and such that, if the initial positions and ve- 
locities belong to this set, the distances of the bodies from each other will remain 
perpetually bounded." 

As well known, such statement solves a fundamental problem considered, for several 
centuries, by astronomers and mathematicians. However, Arnold considered in de- 
tails only the planar three-body case 2 and it appears that his indication for extending 
the result to the general case contains a flaw; compare end of § 1.2, p. 1524 in [7]. 

A complete general proof of Arnold's statement was given only in 2004, when 
J. Fejoz, completing the work of M. Herman, proved the following 3 

Theorem 1 (Arnold, Herman, Fejoz [7, §1.2, p. 1523, Theoreme 1]). 

Si le maximum e = max{mj/mo}j=i n des masses des planetes rapportees a la 
masse du soleil est suffisamment petit, les equations (1) admettent, dans Vespace des 
phases au voisinage des mouvements kepleriens circulaires et coplanaires, un ensem- 
ble de mesure de Lebesgue strictement positive de conditions initiales conduisant a 
des mouvements quasiperiodiques. 

The beautiful proof of this result given in [7] (see also [8]) relies, on one side, on 
the elegant C°° KAM theory worked out by Herman (§ 2-^5 in [7]), and, on the 
other side, on the analytical celestial mechanics worked out, especially, by Poincare 
and clarified and further investigated in Paris in the late 1980's by A. Chenciner and 
J. Laskar in the Bureau des Longitudes 4 and later by Herman himself. 

The invariant tori associated to the motions provided by Theorem 1 , in view of the 
just mentioned KAM tools, are C°°. Now, since the many-body problem is formu- 
lated in terms of real-analytic functions, it appears somewhat more natural to seek 
for real-analytic invariant manifolds. This is the problem addressed in this paper. In 
particular, we shall give a new proof of Arnold's statement, proving the following 

Theorem 2. If e = max{mj/mo}_;=i,...,n is small enough, there exists a strictly 
positive measure set of initial conditions for the (1 + n)— planetary problem (1), 
whose time evolutions lie on real-analytic Lagrangian tori in the Qn-dimensional 
phase space 

n n 

M := {(q,p) € R 6 ^ : q 3 + q k ,Vj^k and = = , 

3=0 3=0 

endowed with the restriction of the standard symplectic form Ylj=o dQj A dpj = 
X^o<j<n dqj k A dpj u. 

l<fe<3 

'The integer n in Arnold's statement corresponds to the above (1 + n). 

2 Few lines after the above reported statement in [1, Chapter III, p. 125], Arnold says: "We shall consider 
only the plane three-body problem in detail. [■ ■ • ] In the final section a brief indication is given of the way 
in which the fundamental theorem of Chapter IV is applied in the investigation of the planetary motions in 
the plane and spatial many-body problems." . 

3 For a more detailed statement, see, footnote 6 below. 

4 Compare, e.g., the Notes Scientihques et techniques du Bureau des Longitudes S 026 and S 028 by, 
respectively, Chenciner & Laskar and Chenciner. 
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Remark 1. Let us collect here a few observations concerning the above statements 
and respective proofs. 

(i) The proof of Theorem 2 given below is similar in strategy to that in [7] but 
technically different and it is based on an analytic (rather than smooth) KAM 
theory for properly-degenerate Hamiltonian systems (see, also, point (iv) be- 
low). On the other hand, it is conceivable - notwithstanding the presence of 
strong degeneracies (see point (ii) and (iii) below) - to prove regularity and 
uniqueness results for the planetary problem so as to deduce that the invariant 
tori in [7] are indeed analytic, a fact which does not follow from our proofs 5 . 

(ii) The evolution equations (1) are Hamiltonian and admit seven integrals, namely, 
the Hamiltonian (energy) H := £"=o fejj - £o<*<i<n the three 
components of the total linear momentum M := £™ =0 Pj and the three com- 
ponents of the total angular momentum C := Yl]=o Pi x Qj> wnere ' x ' denotes 
the usual skew vector product in IR 3 . As a reflection of the invariance of New- 
ton's equation (1) under changes of inertial reference frames, the Hamiltonian 
system associated to the (1 + n)-body problem may be studied on the sym- 
plectic, invariant 6n-dimensional manifold M defined above, where, besides 
the total linear momentum, also the coordinates of the barycenter of the sys- 
tem vanish ("reduction of the total linear momentum"). However, the reduced 
(l+n)-body Hamiltonian still admits, besides the energy, three integrals given 
by the components of C = (C x ,C y ,C z ). Incidentally, such integrals are not 
commuting since, if {•, ■} denotes the natural Poisson bracket on M, one has 
the cyclical relations {C x ,C y } = C z , {C y ,C z } = C x and {C Z ,C X } = C y \ 
but for example |C| 2 and C z are two independent, commuting integrals. 

(iii) The reasons why, notwithstanding the development of KAM theory in the early 
1960's, it took so long to give a complete proof of Arnold's statement are tech- 
nical in nature and are related to the strong degeneracies of the planetary prob- 
lem (degeneracies, which are related to the abundance of integrals mentioned 
in the preceding item). The planetary (1 + n)-body problem is perturbative, 
the unperturbed limit being obtained by considering n decoupled two-body 
problems formed by the Sun and the j th planet. Now, the two-body prob- 
lem in space is a three-degrees-of-freedom problem, but, once it is put into 
(Delaunay) action-angle variables, it depends only on one action (the action L 
proportional to the square root of the semi-major axis of the Keplerian ellipse 
on which the two bodies revolve). Systems of this kind are called properly 
degenerate and standard KAM theory does not apply. This difficulty, however, 
was overcome by Arnold - essentially by refined normal form theory - in the 
case of the planar three-body case, to which he could apply his "fundamen- 
tal theorem" [1, chapter IV]. Indeed, Arnold's approach, in view of Jacobi's 
reduction of the nodes, could be extended [15] to the spatial three-body case 
(n = 2) but not to the general case (spatial, n > 2). Indeed, when n > 2, in 
space, a nice symplectic reduction (corresponding to the reduction of the nodes 
for n = 2) is not yet known. 

Furthermore (but not independently), in higher dimension, there appear two 
secular resonances (see Eq. (52) below), which prevent direct application of 



5 Recent interesting progresses in the study of uniqueness of invariant Lagrangian manifolds appeared in 
[4] and, especially, in [6]; however, as far as regularity is concerned, to the best of our knowledge, the only 
complete proven statement is Theorem 4, §4, p. 34 in [17], which covers only the C°° nondegenerate case. 
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any kind of KAM machinery. We mention that the way we overcome, here, 
this last difficulty is slightly different from that used in [7]: roughly speaking, 
in [7] it is introduced a modified Hamiltonian, which is then considered on the 
symplectic submanifold of vertical total angular momentum; here, we consider, 
instead, an extended phase space by adding an extra degree-of-freedom and 
consider on it a modified nondegenerate Hamiltonian. 

(iv) The main technical tool for us is the analytic KAM theory for weakly nonde- 
generate systems worked out by RuBmann in [16]; the main results of RuB- 
mann's theory (in the case of Lagrangian tori) are recalled in § 2.1 (see, also, 
Lemma 8 in § 2.3.3). The extension of this theory to properly degenerate sys- 
tems is explained in § 2.2 and proved in § 2.3 (which constitutes the longest 
and most technical part of the paper). Finally, in § 3, using several results 
reported in [7], the proof of Theorem 2 is given. 

(v) Finally, we mention very briefly a few questions related to the context consid- 
ered here: 

- Describe, in detail, the motions that take place on the Lagrangian tori. Let 
us clarify this point. From the proof of Theorem 2 given below, in view of 
the indirect argument used ([7, Lemma 82, p. 1578]) we cannot conclude 
that the 'true' motion is quasi-periodic; on the other hand, using different 
arguments, Arnold and Fejoz say that the motion, in the general case, is 
quasi-periodic and takes place on (3n — l)-dimensional tori 6 . Moreover, 
in the spatial three-body case (n = 2) the Lagrangian tori are actually 4- 
dimensional (not 5 = 3n — 1) and the number of independent frequencies 
is 4 (compare [15]). 

- Find a 'good' set of analytic symplectic variables for the general spatial 
many-body problem. 

- Give asymptotic (as e — > 0) estimates on the measure of Lagrangian in- 
variant tori. 

- Apply some of the above result to a subsystem of the Solar system (for 
some progress in this direction, see [5]). 



2 Analytic Lagrangian tori for properly degen- 
erate systems 

In this section we first recall a result due to RuBmann concerning analytic perturba- 
tions of weakly nondegenerate Hamiltonian systems (§2.1) and then show how such 
result may be used to give an analytic version of Herman's C°° KAM theorem on 

6 Arnold ([1, p. 127]): "Thus, the Lagrangian motion is conditionally periodic and to the n 'rapid' 
frequencies of the Keplerian motion are added n (in the plane problem) or2n — 1 (in the space problem) 
'slow' frequencies of the secular motions" . 

Fejoz ([7, p. 1566]): THEOREME 60. Pour toute valeurdes masses m , mi,..., m„ > et des demi grands 
axes ai > • • • > a„ > 0, il existe un reel e > tel que, pour tout e tel que < e < e , le Hot de 
l'hamiltonien F (defmi en (28)) possede un ensemble de mesure de Lebesgue strictement positive de tores 
invariants de dimension 3n — 1, de classe C°° , quasiperiodiques et e-proches en topologie C° des tores 
kepleriens de demi grands axes (ai, a„) et d'excentricites et d'inclinaisons relatives nulles; de plus, 
quand tend vers zero la densite des tores invariants au voisinage de ces tores kepleriens tend vers un. 
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properly degenerate systems (i.e., nearly-integrable systems, which when the pertur- 
bation parameter vanishes depend on less action variables than the number of degrees 
of freedom). The statement of the analytic theorem for properly degenerate systems 
is given in § 2.2 and its proof in § 2.3. 



2.1 RiiBmann's theorem for weakly nondegenerate systems 

We start with fixing some notation. 

• If a, b e R" then (a, b) := X)"=i an d M : ~ M2 : ~ ( a ' a ) 5 ' 

• if g is a /i-times continuously differentiable function (p, e N) from an open set 
£> C K" to M m , the /x-th (tensor) derivative of g in b 6 _B is denoted by (ai, . . . a M ) — > 
d^g(b)(ai, . . . , a M ), otj 6 R n ,j = if 01 = ••• = a p , we shall write 
dVg\b){ay 

• l<9 M .gO)| :=max aeR n ; | | =1 |a' i 5(6)(a,...,a)|and|^fif| A := swp beA |<9^(6)|; 

• C M (-B, R m ) will denote the Banach space of all ^i-times continuously differentiable 
functions g : B — > W n with bounded derivatives up to order y,, endowed with the 
norm \g\% = sup <„< M \d u g\ B < 00. 

The key notion of nondegeneracy is the following. 

Definition 1 (RuBmann nondegeneracy condition 7 ). A real-analytic function 
oj : y G B C R n — > wfo) = (wi(y), . . . , w m (y)) G M m 

w called R-nondegenerate i/5 « a non-empty open connected set in W 1 and if for 

any c= (c±, . . . , c m ) G M m \ {0} orce has 



(c, cj) := ^ CiUJi / 



i=i 



or equivalently if the range uj(B) ofuj does not lie in any (to— 1)- dimensional linear 
subspace ofW 11 . We call oj R-degenerate if it is not R-nondegenerate. 

The following lemma is a simple consequence of R-nondegeneracy and analyticity: 

Lemma 1. Let u : B C W 1 — > M. m be R-nondegenerate. Then for any non empty 
compact set K. C B there exist numbers p,o = /Uo(w, JC) G Z + and (3 = /3(u, /C) > 
such that 



max 
o<m<mo 



dS(cMv)Y 



>f3 , V cGS m -\ Vy G/C (2) 



where S m ~ l := {c G R m : \c\ 2 = 1}. 



For the proof see Lemma 18.2 on page 185 of [16]. 
In view of Lemma 1 one can give the following 

7 This terminology seems to be, nowadays, standard (see, e.g., [18]); however many authors, besides 
RuBmann, contributed to its formulation; among them: Arnold, Margulis, Pyartli, Parasyuk, Bakhtin, 
Sprindzhuk and others. 
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Definition 2. Let JC and B as in the preceding lemma and let uj : y € B — > W 71 be 
a real-analytic and R-nondegenerate function. We define fio(uJ, JC) £ Z + , the index 
of nondegeneracy of uj with respect to JC, as the smallest positive integer such that 



(3 := min max 

yeK, ceS m - 1 0</i<W) 



> o. 



(3) 



The number j3 = (3(uj, JC) is called amount of nondegeneracy of uj with respect to 
JC. 



Remark 2. If a real-analytic function uj : B —> R m admits the existence of /j,q and 
f3 as in (3), for some compact JC containing an open ball, then it is R-nondegenerate. 

The following result, which concerns the existence of maximal (Lagrangian) tori 
only, is a particular case of the main Theorem in [16], where also lower dimensional 
tori are treated 8 . 

Theorem 3 (RuBmann, 2001). Let y be an open connected set of W l and T n the 

usual n-dimensional torus R n /27rZ ra . Consider a real-analytic Hamiltonian 

H(x,y) = h(y) + P{x,y) 

defined for (x,y) £ T" x ^ endowed with the standard symplectic form dx A dy. 
Let K, be any compact subset of y with positive n-dimensional Lebesgue measure 
meas n JC > and fix < e* < meas n IC. Let A be an open set in C n /27rZ n x C n 
on which H can be analytically extended and such that T n x tC C A. Assume that 
the frequency application uj := V/i is R-nondegenerate on y, m let jjl be any integer 
greater or equal to fj,o(uj,]C) (the index of nondegeneracy of uj with respect to JC) 
and let f3 be as in (3) with p,Q replaced by fi. 

Then, for any fixed r > nfi, there exist €q = eo(e*, n, p,, [3, r, uj, JC) > and 7 = 
7(e*, n, r, p, f3, uj,JC) > such that if 

\P\ A := sup \P\ < e (4) 
A 

the following is true. There exist a compact set 

JC* C JC with meas ra JC* > meas„ JC — e* (5) 

and a Lipschitz mapping 

X : (b, £, r]) g JC* x T n x U — > T n x y , 

where U is an open neighborhood of the origin in W 1 , such that: 
( i) the mapping 

(Z,n)^(x,y) = X(b,Z,r,) 

defines, for every b € JC*, a real-analytic symplectic transformation 9 close to 
the identity on TP x U; 



8 See Theorem 1.7 in [16, p. 127]. We refer to [13] for more details on how to obtain Theorem 3 from 
the general results in [16]. 

9 I.e. it preserves the symplectic form dx A dy. 
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( ii) the map 

(b, G £* x TT — > X (6, := £, 0) 
is a bi—Lipschitz homeomorphism; 
(Hi) the transformed Hamiltonian H* := H o X is in the form 10 : 

H*(b,^ V )=h^b) + (^(b),r ] )+0(\r ] \ 2 ) 

for every b G K* and (f , rj) G T n x 
f ivj the new frequency vector oj* satisfies for all b in JC* the Diophantine inequality 

\(k,^(b))\>^, VfceZNjo}. (6) 

\ K \2 

Remark 3. (i) From Theorem 3 we immediately obtain that for any b G /C* the 
n-dimensional tori 

% := X (b,T n ) 

are invariant for H and the ii-dynamics is analytically conjugate to £ — >■ 
£ + co*(b)t. Furthermore, as it follows from (5) and point (ii), the measure 
of Ub(zfc*Tb is proportional to (meas n /C — e*)(27r) ra and, hence, tends to the 
full measure linearly when e* tends to 0. 

(ii) A (technical) difference between RuBmann's Theorem and the formulation 
given above in Theorem 3 is the choice of p, as any integer greater or equal 
than the actual index of nondegeneracy of u, while in [16] fi is chosen equal 
to the index of nondegeneracy of uj. In fact, it is easy to check 11 that RuB- 
mann's theorem holds in this slightly more general case, which will however 
be important in our applications. 

Another difference of Theorem 3 above with respect to RuBmann's original for- 
mulation, concerns the way the small divisors are controlled. RuBmann uses a 
very general approach based upon "approximation functions"; however, such 
approach is too general for our application and cannot be applied directly. Nev- 
ertheless, it is easy to follow a more classical approach 12 based upon Diophan- 
tine inequalities of the form (6), which will be good enough for the application 
to properly degenerate systems; compare also remark 5, (ii) below. 



2.2 A KAM Theorem for properly-degenerate systems 

Let d and p be positive integers; let B an open set in M. d , U some open neighborhood 
of the origin in M 2p and e a "small" real parameter. Consider a Hamiltonian function 
H e of the form 

H £ ((p, I, u, v) = h(I) + ef(ip, I, u, v) , (7) 

real-analytic for 

(if, I, (u,v)) £T d xBxU=:M 

10 Here and in what follows f(rj) = 0(g(rj)) means that there exists a constant C such that |/(?7)| < 
C\g(rj) \ for small enough rj. 
"See, e.g., [13]. 
12 Compare, again, [13]. 
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where M is endowed with the standard symplectic form 

dip A dl + du A dv . 
The "perturbation" / is assumed to have the form 

f{<P,I,u,v) = f (I,u,v) + fi((p,I,u,v) , f Jd fi(ip,I,u,v)dip = 

fo(I,u,v) = foo(I) + ^%(/)^ T ^ + O (\(u,v)\ 3 ) . 

3=1 

Observe that the Hamiltonian h + e/o possesses for every / G B the invariant 
isotropic (non-Lagrangian) torus 

Tf ■= T d x {/} x {0} C M 

with corresponding quasi-periodic flow 

<p(t)=(d I h(I) + ed I f O0 (lj)t + <p I(t)=I («(*),«(*)) =0. 

The purpose is to find Lagrangian invariant tori for H e close to [d + p)-tori of the 
form 

T^J = T d x {7} x {(«, «) € M 2p , \(uj,Vj)\ 2 = 2 Wj , V j = 1, . . . ,p} (9) 
for I in B and w € (M + ) p small. 

Theorem 4. Consider a real— analytic Hamiltonian function H e as in (7) and (8), 
and assume that the "frequency map " 

I £ B — > (u(I),Sl(I)) := (V/i(J), fii (/),... ,fipC0) £» d xR p (10) 

« R-nondegenerate. Then, if t is sufficiently small, there exists a positive measure set 
of phase space points belonging to real-analytical, Lagrangian, H e — invariant tori, 
which are close to Tj^ p as in (9) with Wj = 0(e); furthermore, the H e -flow on such 
tori is quasi-periodic with Diophantine frequencies. 

Remark 4. (i) The above Theorem may be viewed as the real-analytic version 
for Lagrangian tori of the C°° KAM Theorem by M. Herman contained in [7] 
(see, in particular, Theorem 57, page 1559). Under stronger nondegeneracy 
assumptions the above theorem corresponds to the "Fundamental Theorem" in 
[1]. 

(ii) The word "properly-degenerate" refers to the fact that for e = the Hamilto- 
nian Hq depends on d action variables, while the number of degrees of freedom 
is d+p > d. In particular, the tori constructed in Theorem 4, as e — > 0, degen- 
erate into lower dimensional (non Lagrangian) tori Tf. 

(iii) The natural symplectic variables for the KAM theory of the Hamiltonian H € 
are (I, ip) and (rather than the cartesian variables (q,p)) the symplectic action- 

angle variables (w,Q, where Wj = 3 2 3 for j = 1, . . . ,p and Q is the angle 
of the circle 13 \wj\ = const. Indeed, Theorem 4 has, in terms of such variables, 

Compare Eq. (27) below, where w is related to p by w = p° + p. 
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a natural reformulation, which gives a deeper insight into the structure of the 
invariant tori 14 : 

Theorem 5. Let H e (<p, I, £, w) = h(I) + ef((p, I, (, w) be real-analytic for 
(<p,I,£,w) G T d xBxJPx{w G W : < |itfj| < r} =: Mfor some open set 
B Cl' i and < r; M. is endowed with the symplectic form dipAdl+d^Adw. 
The perturbation f is of the form f = fo(I, (, w)+fi with f± having vanishing 
(p-mean value over T d ; furthermore /q has the form /o = foo(I) + (£l(I), w) + 
o(\w\). Then, if the frequency map I G B — > (u, Q) := (djh(I),Q(I)) G 
R d x MP is R-nondegenerate, and e is small enough, there exists a positive 
measure set of phase space points belonging to real— analytical, Lagrangian, 
H t -invariant tori, which have the following parametrization: 

' <p = 6 + <p(9,?P) 
1 = 1 + 1(0,4) 

< 

, w = w + w(9, ip) 

where w is a constant vector of norm 2e and <p,I,£ and w real— analytic func- 
tions for (9,i/>) G T d x Tf (with range, respectively, in T d ,R d ,JP andW) 
with 

' /=o( e (i oge -r (To+1) ) 

< w = 0(e^) 
><P,C = 0(e) 

for suitable v > 4 and 15 To > d+p. Moreover, if(u>, f2) is the above frequency 
map, the H £ -flow on such invariant tori is conjugated to 

(0,ip) — ► (9 + u>t,ip + ent S ) 

for a suitable Diophantine vector (Q, (l) satisfying 

\u-u)\, \n-n\ = 0(e) . 



2.3 Proof of Theorem 4 

First of all, let us introduce some notation and make quantitative the assumptions of 
Theorem 4. 



14 In reformulating Theorem 4 in terms of the variables (£, w) we shall often use the same symbols 
used above. The Proof of Theorem 5 will not be explicitly given since it follows easily from the proof of 
Theorem 4. 

15 See equations (16) and (21) below. 
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For S > 0, d G N, A C R d or C d we denote 



^(xo,(5) 
D d (x ,6) 

A + 5 



= {xeR d 
= {xeC d 
= {xeC d 



\x-x \<5}, (x G R rf ) , (11) 
\x - x \ < 6} , (x G C d ) , 
\Imxj\ < 5, Rexj G T, V j = 1 . . . d} 
[jD d (x,5) 

• We may assume that H e in (7) and (8) can be holomorphically extended for 

I, («, v)) eT d a x(B + r ) x (W + n) =: . (12) 

In particular H e is real-analytic on T d x B d (Io, s) x B 2p (0, ri) for any Iq in 
£> and s < tq. Moreover, we shall denote 

M :=V( sup \f k (I,u,v)\)e^ (13) 
k&d V(B+r )x(W+r 1 ) / 

as the "sup-Fourier" norm of / and let 

Mi:= sup |(o;(J),0(J))| . (14) 

The proof of Theorem 4 is based upon two preliminary steps: 

1 computation of a suitable normal form for H e ; 

2 quantitative estimates on the amount of the nondegeneracy of the normal form. 

2.3.1 Step 1: Normal forms for properly-degenerate systems 

Proposition 1. Fix an integer v > 4. Then, there exists m > d (depending on cu), 
and, for e small enough, a point Iq G B and a real-analytic canonical transforma- 
tion 16 <J? e such that the following holds. Let 

S :=o((lo ge - 1 )- m ) (15) 

then B d (Io, s) C B and <3? e : r, £, p) — > (92, /, u, v) satisfies 

$ e : T d x B d (0, s/5) xfx £ p (0, e) — ► T d x £ d (/ , s) x W 
and H e := H e o Q e takes the form 

H e (#,r,(,p) = N e (r,p; P °) + ^P e ^,r,C,p;p ) (16) 

with 

N e := h(I + er)+g(I + er) + ^n(I + er)-(p + ep) + Q eyIo+er {p° + ep) (17) 

a«<f: p° in « some /?o/«? having euclidean norm 2e; Q ei / + er is a polynomial 

of degree v — 1 starting with cubic terms; g , f2 araci P e are real-analytic functions. 
Furthermore, one has 

sup |6(J + er) - n(/ + er)| = O (e(log e- 1 ) 2 ™ -1 ) . 



^Symplectic up to rescalings. 
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Proof of Proposition 1. We start by recalling a measure theoretical result due to 
Pyartli (see [14] or [16, Theorem 17.1]): 

Lemma 2 (Pyartli). Let JC C M d be a compact set; let 9 G (0, 1) define ICg := 
U ye icB d (y, 9). Let g : ICg — > R be a real-analytic function satisfying 

mm max \d u g(y)\ > p 

y&K, 0<f</io 

for some ft > 0. 77iera ?/iere ex/^?^ C = C(p,o, (3, d, /C, 0) smc/i 

meas d {y G /C : | 5 (y)| < i} < C|<7|£° +1 i^ 
/or any < t < ^f^- 
Pyartli 's Lemma implies: 

Lemma 3. Let K. be a compact set with positive d-dimensional Lebesgue measure 
and let let < e* < meas^ K; let u : ICg — > R. d be R-nondegenerate and let p,Q 
and f3 be its index and amount of nondegeneracy with respect to IC. Let us denote 
by D d Q T0 the set of Diophantine vectors in R rf with Diophantine constants 70, To, i.e. 
the set 



V. 



70, TO 



<ER d : \(oj,k)\ > 



]^». V*GZx{0} 



Then, if is sufficiently small and To > dfio one has 

measrf (jCC\V d ^ > meas^/C-e* 



(18) 



Proof of Lemma 3. First of all observe that for any m G Z+, a G R d , k G 7L d \ {0} 

and 6 G /Qi 

^Mfe),*!*^)^" 1 ) = (^"M&Xa™),*!^ 1 ) < |a m w(6)(a m )| ; 
taking the sup over \a\ = 1, 6 G Kg and < m < /j, we get 



< \w\fc e < 00 



for any G Z+. Now we use this last inequality and Theorem 2, assuming 70 < 



2(10+2 



, to estimate 



meas, (k \ P^) = meas, |J {& G K : |<o;(6),fc)| < | 



fcGZ d \{0} 



< ^ measrf < b G /C : 
fcgZ d \{0} ^ 



< 



7o 



|*| 



TO + l 



< C( / uo,/3, ( i,/C,0)|^|^ +1 7 o" () E 



TQ + l 

fcez«K{o} |*| M0 
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Since tq > djiQ this last sum converges and one has 

meas, (/C \ < ;T0 ) < 



for a suitable C = C(/Uo, (3, d, /C, 0, oj, To). Choosing 70 < (C 1 e*) Mo we obtain 
estimate (18). ■ 

Now consider the real-analytic Hamiltonian H e in (7) and (8). Let v\, v 2 > 4 be two 
integers to be later determined and set 

^1 := -(^i-l) log -J- (19) 

where Mq is defined by (13). Lemma 3 and the R-nondegeneracy of u assure the 
existence of Iq G B such that ui(Iq) belongs to ^ 0;ro (for suitable 70 and To). Then, 
from Taylor's formula it follows that 17 

|w(J) • k\ > qi > , V fc € Z d , < \k\ l < K x , V / € £> d (/ , s) (20) 

with 18 

s-O^loge- 1 ) - ^) and ai := O ((loge^P') . (21) 
Furthermore, we can assume that there exists a 2 (independent of e) such that 

-k\>a 2 >0, V k G IT, < \k\ x < v 2 , V / e L> d (/ , s) . (22) 

Next, we want to average ff e over the "fast angles" ip up to order Z7. To do this we 
shall apply the following classical "averaging lemma", whose proof can be found in 
[3, Appendix A, p. 110]. 

Lemma 4 (Averaging Lemma 19 ). Let H e , Mq, a, a\ and s be as above. Assume 
(20) holds with K\ as in (19). Then, if e is small enough, there exists a real-analytic 
symplectic transformation <J>* : ((p, I, u, v) — > (p, I, u, v) mapping 

Mi := T| x D d (I , J) x D 2p (0, % M := x D d (I , s) x D 2p (0, n) 
that casts H t into the Hamiltonian 

R\ := H e o $1 = h + e/ + 9 + f 
where g = g(I, u, v) and f satisfy 

sup \g\ < , sup l/l < (eM r (23) 

Mi sa l Mi 

for a suitable C = C(a, t>\). 



17 ForA: = (k\, k 2 , ■ ■ ■ , kd) € Z d we denote \k\ x := Y^ =1 | fc^ | ; recall also the definition of complex balls 
D d m(\\). 

18 This means that we can take m = tq + 1 in (15). 

"Lemma 4 can be immediately derived from Proposition A. 1 in [3] with the following correspondences: 
ct\ = a for ct\ as in (20), K\~K for K\ as in (19), eM = e for M as in (13), s — r, d for s as in (15) 
and (20), {0} = A and ef(I, <p, u, v) in (7) is just f(u, ip) in [3]; as a result one has that e/o + g and / are 
respectively given by g and /* in [3] with estimates (23) holding in view of the previous correspondences. 
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Thus, if we set g =: e 2 g and / =: e Ul f, using (8) and (23) we have 

He{<pJ,u,v) = h(I) + e f (I,u,v) + eg(I,u,v) + e Ul f(ip, I,u,v) 

f (i,u,v) = f 00 (i) + j: p j=1 n j (i)^ + o(\uM 3 -J) ■ 



(24) 



From equation (24) we see that the application of averaging theory may cause, in 
general, a shift of order e of the elliptic equilibrium, which, before, was in the origin 
of IR 2p . Therefore, we focus our attention on the Hamiltonian function /o + eg with 
the aim to find a real-analytic symplectic transformation restoring the equilibrium 
in the origin. An application of the standard Implicit Function Theorem yields the 
following 20 : 

Lemma 5. Let M2 ■= T| x D d (I , s/4) x D 2p (0, n/4); then, provided e is suffi- 

7 

ciently small, there exists a (close to the identity) real-analytic symplectic transfor- 
mation 

$ 2 : (x,y,p,q) G M 2 — ► (<p,I,u,v) G Mi 
such that H 2 := Hj o <fr 2 is of the form 

H 2 {x,y,p,q) = h(y) + eg(y,p,q) + e Ul f(x,y,p,q) 

with d p g(y, 0, 0) = = d q g(y, 0, 0), g and f real-analytic on M.2- 

Now, we need to control the frequencies associated to the modified Hamiltonian 

$M,0): 

Lemma 6. If e is small enough then the eigenvalues of the Hamiltonian g, i.e. 
the eigenvalues of 21 J2pd 2 p q ^g{y, 0, 0), are given by 2p purely imaginary functions 

±ifl\, . . . , ±ifl p verifying 

sup \n(y) - n(y)\ = O (e(lo ge - 1 ) 2ro+1 ) (25) 

ydD d (h h s/A) V 7 

for ro as in (21). 

Proof of Lemma 6. Consider the quadratic part of g, that is the real-analytic 
2p x 2p symmetric matrix A(y) := <9? \g(y, 0, 0). Using the construction of 3> 2 
in Lemma 5, g = (/o + eg) o estimate (23) together with g = e 2 g and the defi- 
nition of s and «i in (21), equation (24) for /o and Cauchy's estimate for derivatives 
of analytic functions, one has 

A(y) = diag (fi^y), . . . , n p (y), Q^y), %{y)) + O (e^oge- 1 ) 2 ^ 1 ) . 

20 Since ilj ^ for every j in view of (22), we can apply the Implicit function Theorem to obtain, for 
small enough e, the existence of two functions w = u (I, e) and u = v (f e) which are real-analytic 
for / e D d (I , s/4) and such that Vu,v(fo + uo, vo) = 0. Furthermore, using (23) together with 

g = e 2 g and (21), one has u 0l v — (9(e(loge _1 ) TQ+ ). The symplectic transformation in Lemma 5 is 
then generated by x ■ <p + (p + Uo(x, e)) ■ (v — vq(x, e)) . 

21 J 2p denotes the standard 2p x 2p symplectic matrix. 
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Since Clj / $7^ for j ^ k, an application of the Implicit Function Theorem 
tells us that the eigenvalues of g (that a priori might have non-zero real part) are 

O ^e(loge _1 ) 2r ° +1 ^ close to Now, asitis well known, eigenvalues of Hamil- 

tonians always appear in quadruplets ±A, ±A; thus, from the simplicity of the eigen- 
values of g (holding for e small enough) one has that its eigenvalues are purely imag- 
inary as claimed. ■ 

By normal form theory (see corollary 8.7 of [2]) we can find a real-analytic sym- 
plectic transformation 0(e)-close to the identity 

$ 3 e ■ (x,y,M e M 3 := Ti x D d (I Q , S -) x D 2 P(0, ^) — > (x,y,p,q) € M 2 

with y = y and such that the transformed Hamiltonian function Hf := H 2 o 
which is real-analytic on M 3 , has the form 

p 

H*(x, y,p, q) = h{y) + eg {y) + | ^ (p 2 + g?) 

i=i 

+ m(y,P, q) + ( Vl h{x,y,p, g) 

where ^ := 0, 0), f 3 ■= f o ^ and # 3 := g 3 o $ 3 e verifies 

sup |&(&f>,9)l < C|(M)| 3 V (p,q) & D 2 P(0,n/5). 

yeD d (I ,s/5) 

Now let g2(y,P, q) ■= \ Ya=i ^iiPi + <2?)> we want to P ut 92 + e<73 into Birkhoff's 
normal form up to order v^. In view of inequalities (22) and (25), provided e is small 
enough, we have 

My) -k\>Y> V A; € Z p , < \k\ 1 < V2, V y e D d (/ , s/5) . (26) 

By Birkhoff's normal form theory 23 , one obtains easily the following 

Lemma 7. If inequality (26) is satisfied, then there exist < r* < r[ < r±/5 and a 
real-analytic symplectic diffeomorphism : (9, r, u, v) — V (x, q) mapping 

Ma := Tl x D d (I , J) x D 2p (0, r*) 3 M 3 := Ttl x D d (/ , J) x D 2p (0, r[) 

leaving the origin and the quadratic part of invariant, such that (9,r) = (x,y) 
and Hf := H 3 o ^ is of the form 

v 

Hf(e,r,u,v) = h(r) + eg a (r)+ e -Y J ^(r)(u) + v 2 ] ) + 

3=1 

+ eQ+(r,u,v) + eR+(r,u,v) + e Ul f A (6, r, u, v) 

where: 
From (22). 

See, e g., [10, Theorem 11, p. 43] or [13, section 3.4] for a quantitative version. 
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• is a polynomial of degree [^] in the variables I = . . . , I p ) having the 
form 

<fi(r), J) + \{T{r)I, /> + ..- with Ij := \{u) + vf) 

with T(r) a 2p x 2p real-analytic matrix; 

• R*is a real-analytic function verifying \R+(r, u,v)\ < C\(u, v)\ U2+1 for every 
(u, v) € D 2p (0, r*) and r € D d (I , s/5); 

• /4 := /3 o is real-analytic on M4. 

We may conclude the proof of Proposition 1. Following [7, pp. 1561-1562], we pass 
to symplectic polar coordinates in order to move R+ to the perturbation of Hf with 
the help of a rescaling by a factor e. Let p° = (p§, . . . , p®) in be sufficiently 

close to the origin; consider, for a suitable a+ > 0, the real-analytic symplectic 
transformation ^ : (9, r, (, p) — » (9, 1$ + r, z) mapping 

M 5 := Tl x L> d (0, J) x x (0, ^) ^> M 4 
8 5 2 

where 

Zj = Uj + ivj := ^ 2 (//! • p,) , . (27) 

The transformed Hamiltonian function H\ := Hf o real-analytic on M.5, as- 
sumes the form 

p 

H 5 e (9, r, C, p) = h(I + r) + eg {h + + | E ^( J ° + 0(p° + P?) 

i=i 

+ eQ/ 0+r (p° + p) + eR(I + r, C, p° + p) + r, C, P! P°) 

where 

• Qi + r '■= Q* o is a polynomial of degree [^] with respect to p° + p, 
depending also on Iq + r; 

• i? := i?* o $;> verifies 

|i?(/ + r,C,P° + p)| <C|p°|^ 
for every p € L> 2p (0, |p°|/2), r G Z) d (0,s/5) and ( G 1?/, 

• /5 : = /4 is real-analytic on M5. 
Now, let A e be the homothety given by 

A e : (9,r,C,p) —> (9,er,(,ep). 

Even though ^4 e is not a symplectic map it preserves the structure of Hamilton's 
equations if we consider the Hamiltonian function Hf := ^H^ o A e . Explicitly we 
have 

Hl(9,r,(,p) = h(I + er) + go(I + er) + ^n(I + er)-(p° + ep) 
+Q eJo+er (p° + ep) + R(I + er, ep, C; p°) 
+e^- 1 M9,r,(, P ;p°). (28) 
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"2 + 1 



where / 6 := / 5 o A € . Now we fix p° G with |p°| = 2e so that < Ce~ 

Thus, if we choose v\ and so that 



v x -\ = 



v 2 + 1 



:= v 



(29) 



we may write 

R{h + er, ep, C; P°) + e^/e^, r, C, p) =: e"P e (0, r, C, p) (30) 
for a suitable function P e real-analytic on T£ x L> d (0, s/5) x T£ 4 x WO, e) . 

8 

We have proved Proposition 1 with H e = Hf as in (28), (30). ■ 

2.3.2 Step 2: Amounts of nondegeneracy of the normal form 
Proposition 2. Let N e be as in (17). Ife is small enough, the frequency map 



^ t : (r,p) <E B d (0,s/5) x £ p (0,e) 



Or op 



is R—nondegenerate. 



Moreover, let p and j3 denote respectively the index and the amount of nondegeneracy 
of the unperturbed frequency map (10) with respect to a closed ball B d (Io,t) C B, 
for some t > independent ofe. Then, if we define JC e := B d (0, s/W) x B p (0, e/2) 
and let p t denote the index of nondegeneracy of^ t with respect to K. e and 



f3 e := min min max 

ce Sd+p-i (r,p)<EK e 0<M<£ 



^,p)lMe>f 



one has 



p e < p and $ e > 



(31) 



Proof of Proposition 2. From (17), it follows that the frequency map of N e is given 
by 

* e (r, p) = (oo(I + er) + O(e), Ul(h + er) + 0(e 2 )) 

and it is real-analytic on D d (0, s/5) x D p (0, e). Using (25) one has 

* e (r,p) = (w(J + er) + O(e), | (fi(/ + er) + O(e))) . (32) 

Now, let p: G N + and (3 > denote respectively the index and the amount of nonde- 
generacy of ^ := (ou,Q) with respect to B d (Io,t), for some positive t independent 
of e. Set 

*o(r) := (w(/ + er), fi(/ + er)) , (33) 
/Co := .B d (0, s/10) and use definition 2 to get 



mm max_ 

r£/Co 0<fi<fi 



^|<c,*o(r))r 



> e^/3 > 



16 



for every c G 

Next, denote by \£ e the real-analytic function over D d (0,s/5) x D p (0,e) obtained 
multiplying the last p component of ^ e by a factor 2/e. Then, observe that equations 
(32) and (33) imply * e (r, p) = ^ O (r)+O(e). Therefore, denoting /Ci := B p (0, e/2) 
and assuming e small enough, one has 



f3 t := min max 

(r,p)elC xK 1 0<^<p 



^ } Kc,* e (r,p)>| s 



e^/3 



for every c G 1 . 

Now, if we write \I/ e = , Vl/( 2 ) ) el^x IF, from what observed before, it results 

*e=(^ 1} ,|^ 2) ) • 

Define for c = (ci, c 2 ) G M d x IIP with |c| = 1 the function 



/(r, p, ci,c 2 ) 



max 

0<p<£ 



= max 

0<p<£ 



(np) 



furthermore set 



|2 , 6 I |2 

ci| + jM 



and c\ = c\t e 1 , c 2 = ec 2 (2t e ) 1 so that |(ci, c 2 )| = 1. Then one has 



max 

o<p<a 



<9 M 

(r,p) 



= / (r,p,ci,|c 2 ) = i 2 / f r,p, 



ci e c 2 



L'2L 



e 2 e^ +2 /3 
> -£f(r,p,a,ci) > — — > 



and it follows immediately 



mm max 

(r,p)G/C x/CiO<p<Ai 



P P+2 



> 







> 



for every c G 5 d+ * ,_1 . Since /C x /Ci = B d (0, s/10) x B*>(0, e/2) = /C e we have 
verified (31). In view of remark 2 we also conclude that \f e is R-nondegenerate on 
i? rf (0, s/5) x £? p (0, e), provided that e is small enough. 

Proposition 2 is proved. ■ 



2.3.3 Conclusion of the Proof of Theorem 4 

We want to apply RuBmann's Theorem 3 to the properly degenerate case of H € in 
(16). With Propositions 1 and 2 we are in a position to meet the hypothesis of R- 
nondegeneracy of the frequency application required in Theorem 3. However, the 
"degenerate" case of H e requires that the size of its perturbation is of a sufficiently 
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small order in e. From (16) we see that the size of the perturbation of H e is order e u 
where v can be chosen to be arbitrarily big 24 but independent of e. 

Next we provide an explicit expression for the admissible size of the perturbation in 
RiiBmann's Theorem, i.e. eo in (4). 

Lemma 8 (RuBmann). Let H, y, e*, A and r be as in Theorem 3 and let u := Vh 

be R-nondegenerate (as in the hypotheses of Theorem 3). Consider the following 
quantities: 

1. Let K C y be any chosen compact set; let p, be any integer greater than the 
index of nondegeneracy of uj with respect to K, and let [3 be the "amount of 
nondegeneracy" corresponding to p. 

2. Let € (0, 1) be chosen such that 25 T$ x (JC + M) C A and define C\ : = 
M;c+3i?- Let do be the diameter ofK, i.e. do ■= sup xye i C \x — y\. 

n + l 

3. Let Tq > e r such that the following inequality holds 

1 T 

,rr - ,) . (34) 



4. Define 



and set 



5. Finally set 



7,.r 

Then £q in (4) can be taken to be 



To 



C* := 2^ +1 {l1 + (Ci + 1) (35) 



7 := WCT f ^^- (36) 



-(r+n+l) „ 

f „ - *. (37) 



E 1 := 7 T - (r+n+1) ^ 



E2 := Mr- (38) 



eo:=co^-(miii{£7i,^}) 2 (39) 
for a suitable Co = Co(n, fi). 

Remark 5. The above result follows from [16] by considering the case of maximal 
tori only 26 . More precisely: 



24 Recall (29) and the fact that both v\ and v 2 can be arbitrarily fixed at the beginning of the process 
described in section 2.3.1. 
25 Recall definition (12) 

26 Compare in particular the estimates listed on page 171 of [16]; see, also, chapter 2 of [13]. 
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(i) The maximal case corresponds to the easier case p = q = in [16]. Notice, 
however, that it is not sufficient to substitute the values p = q = in RuB- 
mann's estimates (as ,when p = for instance, many terms in [16, p. 171] 
become meaningless) but, rather, one has to go through the most of Theorem 
18.5 in [16] to get the value of C* and 7 in (35) and (36) and through the first 
part of Lemma 13.4 in [16, pp. 158-161] to get the value of to in (37). 

(ii) (Control of small divisors) In [16, section 1.4] RuBmann introduces a so called 
"approximation function" $ in order to control the small divisors. Our choice 
is to take <1>(T) = T~ T with r > np. Comparing [16, section 1.4], one 
sees that such $ does not verify property 3, i.e. T X Q(T) T -^ ) q for any 
A > 0. However, when we consider H = H e , we will see below that one 
has T = T 0i£ = 0(e~ 2 ); then, equation 14.10.10 together with 13.1.4 and 
inequality 14.10.11 in [16] would cause 0{e u ) to be an inadmissible size for 
a perturbation. Nevertheless we claim that the only decay property which is 
actually needed in RuBmann's Theorem 3 is 

lim T A $(T) = for all < A < np 
so that our choice is perfectly suitable. 



Now we are going to analyze what happens to the estimate in Lemma 8 when we 
consider H e as Hamiltonian function. In particular we are going to show that each 
one of the quantities appearing in Lemma 8 can be controlled by constants involving 
initial parameters related only to H e in (7) and (8) times powers of e. 

We point out that in the application of Theorem 3 with H = H € in (16) we have the 
following correspondences 27 : 

T n = T d xT p , x = (8, C) 

y = y e := B d (0, a/5) x &>{0, e) , y = (r, p) 

(40) 

A = A e := Ti x TTt x D d (0, a/5) x £>p(0, e) 

8 

P = e y P t , N = N e . 

Accordingly to 8. 1 we consider the frequency application of the integrable part of H e , 
that is ty € (r, p) as in Proposition 2. We already proved that ^> e is R-nondegenerate 
for (r,p) € B d (0,s/5) x B 2 P(0,e). Now, in view of 8.1 and the correspondences 
in (40) we need to fix a compact set /C = )C e C A e . For our convenience we take 
K € := B d (0,s/10) x Bf(0,e/2) so that the first inequality in (31) allows us to 
consider 28 p = p, as an integer greater than the actual index of nondegeneracy of ^ e 
with respect to /C e . Also, in view of the second inequality in (31), we can take 

? = (41) 



See Theorem 3, Proposition 1, (12) and (21) for notations. 
Recall Proposition 2 for the definition of p. 
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in (36) and (38). 

Next, we choose i? = i? e := e/16 so that, for e sufficiently small and in view of (15), 
one has K e + 4$ e C .4 e as required in 8.2. Accordingly to Theorem 3 we also need 
to fix a positive number e* < meas^+p /C e . In view of our definition of /C e and (15) 
a suitable choice is given by 

e* = e p+1 (42) 

for e small enough. 

Now, observe that the quantities C\ and do in Lemma 8, point 2 do not cause any 
change in the order in e of the size of the admissible perturbation. In fact, using 
equation (32) and taking e sufficiently small, we have 

Ci = Ci, e := |tf e |jc e +30 e < 



< sup |w(/ + er)| + e sup |fi(J + er)| + 0(e) < 

rGD d (0,s/5) rGD d (0,s/5) 

< sup \u(r)\ + e sup |0(r)| + O(e) < M 1 

reD d (I ih s/5) reD d (I Q ,s/5) 

where Mi is defined in (14). Since the estimate for eo is decreasing with respect 
to C\, we can substitute C\ in (35) and (37) with M\. The estimate for eo is also 
decreasing with respect to do so that when we consider K. = K e we may simply 
replace do by 1. 

Let us now analyze the quantities in 8.3 and 8.4. First of all observe that in view 
of (31) and n = d + p we can fix a priopri an exponent r > (d + p)p, satisfying 
the requirement in 8.3. Furthermore, given the previous choice of # e , inequality (34) 
becomes 

d+p+l 

which can be easily fulfilled, together with To > e t , by choosing 

T = T , e := 1 (43) 

for e sufficiently small. For what concerns the quantities defined in 8.4 we see that 
since # = $ £ := e/16 and the estimate for eo is decreasing in C*, we can choose 

C* = C* := 2 5 ^ +1 V + lf +2 {M x + l)e~ (/m) (44) 

having also used C\ = C\ je < M\. From the fact that we can replace do by 1 
together with equations (41), (42) and (44), one has 

7 = 7e := C1 (M X + l)-f e (A+D 2 +^^ e *f (45) 

for a suitable constant c\ < 1 depending only on |U. Moreover, given once again the 
previous choice of $ = i? e together with equation (43) and the above definition of 
7 e , we can replace to in (37) by 

^ = t ,e = Ci(M! + e (M+l) 2 + ^+2(r+d +P )+3 (4f)) 
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for ci < 1 as above. 



From (38) we see that E\ and E2 have simple polynomial dependence on the quan- 
tities 7, Tq 1 , (3 and to- Our previous analysis shows that when we consider H e as 
Hamiltonian function, these quantities can be replaced respectively by 29 

7, = O ( e (^) 2 +^) , T ;i = 0( e 2 ) , ^ = 0(e) 

f3 e = (e^+ 2 ) , t , e = O (>+l) 2 + ^+2(T + d+ P )+3^ 

Therefore, in view of (39) the size of the perturbation allowed by RuBmann's Theo- 
rem when we consider H = H e , is order e u ° with 30 

vq := 2f + (p + 5)/i 2 + [14 + 4(r + d + p)]fi + 13 . (47) 

In particular we have a condition of the form eo < ce u ° where c is some positive 
constant independent of e and depending only on quantities related to the initial 
Hamiltonian H e , namely p,, (3 and K as in Proposition 2, the Diophantine constant 
r > id + p)p, and M\ as in (14) with tq as in (12). By Proposition 1 we know that 
we can assume the size of the perturbation of H e to be order e u for any fixed integer 
v > 4 independent of e. Thus, by simply taking 31 v > vq, we can apply RuBmann's 
Theorem to H e and obtain Theorem 4 as a consequence. ■ 

3 Proof of Theorem 2 

As it follows from the analysis described in [7, Sect. 6, pp. 1563-1569], the motions 
of (n + 1) bodies (point masses) interacting only through gravitational attraction, re- 
stricted to the invariant symplectic submanifold of vanishing total linear momentum, 
are governed by the real-analytic Hamiltonian 

F = H°(A) + e(i? 1 (A, £, r,, q,p) + H 2 (X, A, £, r,, q,p)) (48) 

where: 

(i) (A, A, f , n, q, p) € T n x (0, oo) n x R n x W l x W l x M n are standard symplectic 
coordinates; 

(ii) Aj = [ijy/Mj a,j, where > are the semi major-axis of the "instantaneous" 
Keplerian ellipse formed by the "Sun" (major body) and the j th "planet", while 

= 1 , Mj := mo + errij , 

efj,j too errij 

too and errij being, respectively the mass of the Sun and the mass of the j th - 
planet; 

'See (45), (43), (41), (46) and recall i? = tf £ := e/16. 

'Using (47) we are able to define the values of v\ in (19) and (20) and V2 in (22) through equation (29). 
Notice that this can be done since v only depends on d, p, r and p,. 
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(iii) the phase space M is the open subset of T n x (0, oo)" x l n x I" x l n x W l 
subject to the collisionless constrain 

< a n < a ra _i < ■■■ < a\ 

n 

and endowed with the standard symplectic form ^ dXj A dAj + d£j A drjj + 

j'=i 

dgj A dpj', 

(iv) := -Fkcp is the Keplerian integrable limit given by 

# - ^ - 2. — 2Af- . 

j = l 3 

describing n decoupled two-body systems formed by the Sun and the j th 
planet; 

(v) the "secular" Hamiltonian H 1 has the form 32 

£2, 2 n 2, 2 

tf 1 = Co + ^aj^p. + + 0(4) (49) 

j"=i j"=i 

where Co, cr.; and depend on A; "0(A)" denotes terms of order greater than 
or equal to four in (£, 77, (and depending on A); 

(vi) H 2 has vanishing average over A G T n ; depend also (in a regular and non 
influential way) on e. 

Remark 6. The variables (A, A, £, 77, 5, p) are obtained from standard Poincare vari- 
ables after a rotation in (£,rj,q,p) needed to diagonalize the quadratic part of the 
secular Hamiltonian; the "eigenvalues" <jj and qj are the the first Birkhoff invariants 
of the secular Hamiltonian; compare [7, pp. 1568, 1569]. 

The frequency map of the planetary Hamiltonian F is given by 

{v\, . . . , u n , (7i, . . . , a n , . . . , q n } 
where the Uj 's are the Keplerian frequencies 



dA i a] A * 



(50) 



It is customary to consider the frequency map as a function of the semi-major axes a 
(rather than of the actions A); we shall therefore call the "planetary frequency map" 
the application 33 



a : a G A 1 — > -j>i, ...,u n ,ai,.. . ,(T n ,Si, •••,?«}€ M 3n (51) 



32 There is a difference of a factor | with the notations used in Ref. [7]. The computations are performed 
in [11]. 

33 Obviously, the property of being R-nondegenerate can be equivalently discussed in terms of the A's or 
in terms of the a's. 
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where 



A := {(ai,a 2 , . . . ,a n ) G R™ : < a n < a n _i < ■ ■ ■ < ai}. 



Clearly, the idea is to apply Theorem 4 to the real-analytic Hamiltonian F in (48) 
with d = n and p = 2n: ((p, I) corresponding to (A, A) here and u corresponding to 
(£, q) and v to (77, p). However, it turns out that the main hypothesis of Theorem 4 
does not hold, namely, the planetary frequency map a is R-degenerate: in fact (up 
to rearranging the (g,p)-variables) one has 








(*j+9) = 0. 



(52) 



The first relation is related to the rotation invariance of the system; the second relation 
seems to have been noticed (at least in this generality) for the first time by Michael 
Herman and is therefore normally referred to as the "Herman resonance". 

The two resonances in (52) are, however, the only linear relations identically satis- 
fied; in fact in [7, Proposition 78, p. 1575] it is proved the following 

Proposition 3. For all n > 2 there exists an open and dense set with full Lebesgue 
measure U C A, where ctj 7^ cti whenever j ^ i and the following property holds: 
for any open and simply connected set V C U, the ctj define 3n holomorphic func- 
tions and if 



for some c 1 6 



a ■ (c 1 , 
, than 

either 



c 2 ,c 3 



v ■ c 



a • c 2 + ? • c 3 







or 



c 1 = , c 2 = , c 3 = (0,..,0,1) , 
c 1 = , c 2 = = c 3 . 



(53) 



In order to remove the secular resonances (52), we consider the following "extended 
Hamiltonian" on 

M := M x T x R 
adding a pair of conjugate symplectic variables 34 (0 p ,p) G T x R: 

2 n 1 

F := F + P - + ep 2 C z with C z := £ (a, - -(£ 2 + r/ 2 + q) +p 2 )) . (54) 

j'=i 

Let us make a few comments. 

(vii) C z is the vertical component of the total angular momentum in Poincare vari- 
ables (compare [12] and also formula (44) in [7]); the form of C z is unchanged 
in the above variables (£, 77, q,p), which are obtained from the Poincare vari- 
ables by an orthogonal transformation. 



I.e., M is endowed with the symplectic form ^ [d\j A dkj + d£j A dr/j + dqj A dpjj + d9 p A 
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(viii) Since C z is an integral for F (i.e., Poisson commutes with F), F and F Poisson 
commutes: 

{F,F}~={F,F} = 

where {•, -}~and {•, •} denote, respectively, the Poisson bracket on M and on 
M; clearly, since F does not depend explicitly upon the angle 6 P , also p is an 
integral for F (and for F). 

This fact will be important later since from Lagrangian intersection theory 
it follows that two commuting Hamiltonians have, in general, the same La- 
grangian tori (see item (x) below for the precise statement). 

(ix) The extended Hamiltonian F may be rewritten as 

F = H° + e(# 1 + H 2 ) 

with 

H° := FKep(A) + ^ , 

n 

H 1 :=C (A)+ / 9 2 ^A i + (55) 
3=1 

£2, 2 n v 2 + a 2 

+ E - P 2 )^ + E (* - + 0(4) . 

3=1 3=1 

Thus, the "slow" action variables 35 are / = (p, Ai, A n ) and the (extended) 
planetary frequency map is given by 

a : (p,a) e AxR* — ► a(p, a) := ((p, v), a, € R n+1 x R n x R n 

with 

oh,- := oj - p 2 , qj := q - p 2 . 

Proposition 3, implies immediately that a is R-nondegenerate: suppose, in 
fact, that 

&■ (foc 1 ),^,^) = pc + v-c 1 + a-c 2 + q-c 3 = 

for some eel and c % <G R n ; such expression is a second order polynomial in 
p and in order to vanish identically have to vanish its coefficients, i.e., 

v ■ c 1 + a ■ c 2 + q • c 3 = , (56) 
c = , 

n 

-E c I + c l = °- ( 57 ) 

3=1 

But then, by Proposition 3 (and because of (56)), one must have one of the 
alternatives listed in (53), which are incompatible with (57). 
Thus a is R-nondegenerate as claimed and Theorem 4 can be applied to the 
extended Hamiltonian 36 F, yielding, for e small enough, a positive measure 



^Corresponding in Theorem 4 to I = Id), d = n + 1; compare also footnote 36 below. 

36 The correspondence with the notation of Theorem 4 being: d = n + 1, p = 2n, H e = F, f = H 1 + H 2 , 

I = (A, p), («, v) = ((£, q), (V,P)), h(I) = H°, foo = C Q + P 2 £" K " = P)> « = 0- 
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set of real-analytic (3n + l)-dimensional Lagrangian tori in M invariant for 
F and carrying quasi-periodic motion with Diophantine frequencies. 

The fact that F is independent of 9 P and that := d p F = (p + AepC z ) 
is constant along F-trajectories (compare point (viii) above) implies immedi- 
ately that the tori T C M obtained through Theorem 4 have the following 
parametrization 

T := { (Zty, e p ), P , p) : (Y>, 9 P ) G T 3 " x t} (58) 

where Z G M and with F-fiow given by 

( j) t ~(Z^,d p ),e p ,p) = (Z(il> + ut,e p + cu*t),6 p + u*t,p) , 
for a suitable vector oj G M 3n , so that (a;, a;*) forms a Diophantine vector in 

M 3n+1_ 

(x) In [7, Lemma 82, p. 1578] the following statement is proved 

If F and G are two commuting Hamiltonians and if T is a Lagrangian torus 
invariant for F and with a dense F-orbit, then it is also G-invariant. 

Thus, since F and F (viewed as a functions on M) commute, the tori obtained 
in (ix) (on which any F-orbit is dense) are also invariant for the flow on hA 
generated by F. Furthermore, the F-flow in M. leaves both 9 p and p fixed so 
that, for any fixed 9 p G T, the 3n-dimensional torus 

T 9p :={{Zty,e p ),e p ,p) : ip G T 3ra j 

is invariant for F. But this means that such tori are invariant also for the F- 
flow in M., finishing the proof of Theorem 2. ■ 

Remark 7. The strategy followed here is similar to that followed in [7] with a few 
differences: first, in [7] p is treated as a dumb parameter and no extended phase space 
is introduced (but an extra argument is then needed to discuss the nondegeneracy of 
the frequency map with respect to parameters and to discuss the measure of the tori 
obtained); secondly, in [7] there is a restriction to a fixed vertical angular momentum 
submanifold, which is not needed here. 
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